function [Vector]=thermal_deformation
fid1=fopen('153.pch','r');
j=1;
Vector=zeros(0,7);
while feof(fid1)==0
    tline1=fgetl(fid1); 
    matches1=strfind(tline1,'$LOAD FACTOR =   1.0000000E+00');
    if length(matches1)==1   
        tline1=fgetl(fid1); 
        while feof(fid1)==0
            Vector(j,1)=sscanf(tline1(1:16),'%d');   
            Vector(j,2)=sscanf(tline1(20:38),'%f');
            Vector(j,3)=sscanf(tline1(39:56),'%f');
            Vector(j,4)=sscanf(tline1(57:73),'%f');
            tline1=fgetl(fid1); 
            Vector(j,5)=sscanf(tline1(20:38),'%f');
            Vector(j,6)=sscanf(tline1(39:56),'%f');
            Vector(j,7)=sscanf(tline1(57:73),'%f');
            j=j+1;
            tline1=fgetl(fid1); 
        end                   
    else       
    end      
end
fclose(fid1);
end